Simple Model for Energy and Force Characteristics of Metallic Nanocontacts* 
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The quantum size oscillations of the energetic properties and the elongation force of the gold 
slabs and wires, isolated and in a contact with electrodes, are calculated in a free-electron model. A 
simple relation between the Fermi energy and the square-potential-well depth is used and tested for 
low-dimensional systems. It is shown that considering the electron subsystem of a slab (or wire) in a 
contact as open one, the contact acts like a sort of electron pump which sucks or pumps out electrons 
from the sample. The effect of the contact potential difference on the elastic force oscillations is 
considered. The calculated amplitudes of force oscillations are in a qualitative agreement with those 
observed experimentally. 
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I. INTRODUCTION 



Recent experimental investigations of gold point- 
contacts by using the atomic force microscope [l| have 
demonstrated a violation of the Hooke and Ohm laws: 
a deformation of the contact leads to a coherent step- 
wise variation of the conductance and force. Theoret- 
ical interpretation of these observations for a deformed 
nanowire was given in the framework of the simple free- 
electron models 0, IE & IE IE • The energetic picture of 
the nanocontact is more complicated, however, than that 
adopted in 0, IE EL IE El First, in the previous considera- 
tions the contact potential difference |8| is not taken into 
account. Secondly, the metallic contact changes its di- 
mensionality with the variation of sample-cantilever dis- 
tance. At first, its shape is rather similar to a slab, to 
resemble a wire at the moment of rupture [E llfjfl . Thus, 
we have to do with a transformation of the 2D open elec- 
tron systems into the ID one. The purpose of this work 
is to investigate the limiting cases of this transition. We 
study the size dependence of the work function, ioniza- 
tion potential and elastic force of the isolated quantum 
slab and wire that are inserted between two electrodes to 
form a contact. We begin with consideration of the en- 
ergetics of an isolated sample. The results are exploited 
later to determine the contact potential difference (CPD) 
for a sample connected with electrodes. 

In order to trace analytically the effect of dimension- 
ality on the energetic and force properties one needs a 
simple model relating mutually Fermi energy, work func- 
tion and the potential barrier. We consider a rectangular 
sample of the volume V = ax ax L, where L denotes the 
dimension along the x-axis. The inequalities a> L and 
a -C L correspond to the geometry of a slab and a wire, 
respectively. The consideration of these two asymptotic 
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limits which simulate the early and late phases of elon- 
gation, allows us to trace the evolution of the energy and 
force characteristics of the 2D and ID metallic structures. 

We have assumed that the energetics of a finite sys- 
tem of bounded electrons can be described by a square 
potential well of the width L (slab) or a (wire) and the 
depth U < 0, 



-U(E F ) = W(E F ) + E F , 



(1) 



where W is the work function and Ep is the Fermi energy 
of a finite sample. 

Let us remember that the potential well of the depth 
U of metal sphere shows the size oscillations The 
Fermi energy of metal film also has a size oscillations 
[l2|. Due to an expansion of the chemical potential of 
spherical cluster of R radius fi = /j,q + fii/R in powers of 
the inverse radius 0] we paid attention on the inequality 



W = Wo - ^ < Wo. 



(2) 



The criterion W < Wo (subscript zero labels the quan- 
tities for a semi-infinite metal) for one-dimensional sys- 
tems was introduced by us [l4| to check the correct size 
dependence of the U{E F ) (by means E F (L)). In order 
to determine the size oscillations in the potential well 
depth U{E F ) we have employed a simple expression for 
the Wq(Efo) provided by [TE|- The latter is based on the 
concept of the image force action 0, 0| and a sponta- 
neous polarization of metallic plasma, to determine the 
distance, from the classical metal surface at which the 
image force begins to act (see also and references 
therein) , 



W (Efo) = 



B 
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(a.u) 



(3) 



J F<d 



where B is the adjusting parameter and r s o is the average 
distance between electrons in the bulk, B — 0.3721 a.u. 
and r s =3.01 bohr radius for gold. 

Adapting the result for the low-dimensional systems 
we replace Wo{E FO ) by W(E F ) and suppose W — > W = 



2 



4.30 eV M 

as L — > oo 14]. In order to get a qualitative 
agreement of calculated force with the experimental data 
we have to assume that Au is monovalent 000 



where is the number of highest occupied subband, np 
equals tie in Eq.Q with the change k — > hp. 

The total number of electrons in a wire is given by 



II. BASIC RELATIONS FOR AN ISOLATED 
SAMPLE 



For the 2D and ID metallic systems the allowed en- 
ergy levels (th e electron kinetic energies) form a quasi- 
continuum |20| . 
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For an infinitely deep potential well of the width a (the 
shape is a cube) this expression reduces to E nxyz = 
h 2 Tr 2 nly z /(2ma 2 ) with n 2 xyz = n 2 x + n 2 y + n 2 . 

For simplicity, we assume that the wave vector com- 
ponents are the solutions of the transcendental equations 
for square potential well: 



k n M = mr — 2 arcsin(/c n /fc ), 



(4) 



where n — 1,2,3,..., and n = n x for M — L and 
n = n v ,n z for M — a, hk° = y/—2mU. One gets two 
relations of identical form, to determine k n and k Uz . 
The inequalities max{n s ,ii z } 3> maxln^} for slab and 
max{fi 5 ,n 2 } <C max{n x } for wire correspond to highest 
occupied level. 

The density of electron states, D(E), is defined by the 
sum X)n XH , 8{E — E nxys ) over the allowed states. Replac- 
ing the three-dimensional summation in the fc-space by 
integration over k y and k z (or over k x ) and summation 
over n x (or over n y and n z ), we get 



D(E) 



nLh 2 



rip, 



and 
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(5) 



(6) 



for the finite slab (or wire, respectively). 

In the equation J5J for a slab Ue is the integer part of 
number, 
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kL + 2 arcsm(k / k°) 



(7) 



where hk = \j2mE. 

In the equation (|HJ for a wire the plus sign in the limit 
of summation indicates that n y and n z run from 1 to 
the maximum value for which the expression under the 
square root is positive. 

Subsequently, the total number of electrons in a slab 
is given by 



N 



2 n F 
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(8) 



N = 2L 
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(9) 



Here N = TiV, where n — 3/ATrr^ is the electron density 
in the bulk of semi-infinite metal. 

Discussing the effect of dimensionality and charging, 
it is useful to analyze ionization potential which is a 
well defined quantity for an arbitrary size of a sample 
pH |22| . The ionization potential which is defined as a 
work needed to remove an electron from a neutral metal- 
lic sample, can be expressed as 



IP = W+ 2C> 



(10) 



where C is the capacitance of the sample. It should be 
mentioned that Wo and IP are the experimentally mea- 
sured quantities, while W has mainly a methodical sense 
because it can be measured only in the limit of C — > 00. 
An extended thin slab or wire of infinite length has an 
infinite capacitance, C — > 00, therefore IP — > W. Equa- 
tion l|10|) can be interpreted as the effect of charging on 
the work function of the neutral finite sample. It should 
be noted that the size correction in W(L) which is simi- 
lar to that of the spherical cluster (see Eq.lO) competes 
with the term e 2 /2C. Since it is impossible to derive an 
analytical expression for the rectangular sample, in or- 
der to estimate the monotonic size dependence of IP we 
apply the well-known formula for the capacitance C of 
the disk of the width L and of the needle of the length L 

m. 

Solving the set of Eqs.Q and (SJ or ©, by using JJJ 
and adapted Eq. @ , one can determine the Fermi energy 
Ep and subband energies for a given L under condition 
V = constant. Consequently, one gets the work function 
W, as function of the Fermi energy of finite sample. 

The elastic elongation force acting on a finite sample 
is given by F = —dE t /dL, where E t is the total en- 
ergy of the sample. Neglecting the temperature effects in 
the adiabatic approximation, from the virial theorem we 
have E t — —K, which means that the energy of a bound 
electron-ion system is negative El • (The virial and 
stress theorems for electron-ion mixtion of solid metal 
representation was formulated in Ref. 25]). Conse- 
quently, 



F = 



dK 

1l' 



(11) 



In expression (|llfl we take into account the potential en- 
ergy of the electron-ion system. This gives aplus sign in 
front of K (an opposite sign appears in 0. where 
the simplest model of nanowire has been exploited) . The 
plus sign follows from the application of the virial theo- 
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The total kinetic energy of electrons are given by the 
expressions 

^ = |^E / tf + *Z.)**KdK = 
k n x 

2 n F 

and 

2L / 2m , 

(E F + 2E ny + 2E nz ), (13) 
for an isolated slab and wire, respectively. 

III. THE EFFECT OF A " POINT "CONTACT 

When a sample is inserted in a contact between the 
reservoirs then its electron subsystem has to be consid- 
ered as an open one, under the condition W{L) — Wo, 
where W and Wo are the work functions for an isolated 
sample and reservoirs, respectively. Due to the contact 
potential difference, S<p, the electroneutrality of the sam- 
ple breaks down and SN electrons from the electron sub- 
system are transferred into reservoirs. In order to de- 
termined 8cj> one can imagine a simple energetic cycle in 
which electronic charges are transferred from a sample 
to infinity and then into electrodes. By expressing the 
ionization potential, IP, of the sample charged by +eSN 
in a form similar to that for the charged spherical metal 
clusters [2(| we have, 



the two faces of a sample the surplus positive charges 
have a similar surface distribution both in the case of 
real ionization of an isolated sample and for a sample 
in contact. (This situation is similar to that encoun- 
tered in single-electron devices 28] ). Then, assuming 
that C c = e8N/5<p, SN < N and A ->■ 0, Eq.fT^ one 
gets 

6</> = (W -W)/e>0. (16) 

For example, the energy spectrum Eni^ of the remaining 
electrons in the slab, N% = N — SN, can be found by 
solving @ for a square potential well of depth 

Ui=U- eS(f>, (17) 

where U corresponds to the isolated slab. Comparing 
(|17|l with and using (|16|) we find that in equilib- 
rium between a sample and the electrodes, the magni- 
tude of Epi is equal to that of an isolated sample (Ep). 
The value of SN can be calculated using Eq.© with the 
change N — > N% and k Ux — > fc„i x . The total kinetic 
energy, K\, of the remaining electrons is defined as ear- 
lier fEa. 1)120) with the changed energy spectrum and the 
number of electrons. The same approach was used for 
the nanowire. In the case of an open system the elastic 
force is determined by a surplus pressure, relative to the 
reservoirs, multiplied by the area of the contact: 

dfl , . 

where the size-dependent part of the grand potential 
is given by 

Q = E-l + WqN!, (19) 

and Ei = —K\. 



IP — En-SN-A — Ejy^SN — 

WA+^-[(SN + A) 2 -SN 2 }, (14) 

where — eA is the part of the electronic charge which 
leaves the charged sample. The electron affinity of this 
electronic charge, transferred to reservoirs is EA — Wo A. 
The equilibrium condition IP — EA = leads to 

e 2 

W - W-— -(2SN + A) = 0. (15) 

It should be noted that the magnitude of A can be in- 
finitesimally small (non-integer), because in a contact 
(open system), the residual electron can be transferred 
only partially (i.e., there is a finite probability that it 
could be found both in a sample and in a reservoir). 
Thus, we can ascribe to SN a continuous value. We 
also suppose that C of rectangular sample appearing 
in Ea. l|14|) corresponds to the total capacitance C c of 
both contacts. This is justified by the fact that near 



IV. RESULTS AND DISCUSSION 

The calculations were performed for the set of isolated 
Au slabs and wires and then for the ones in contact with 
reservoirs. It allowed us to calculate the contact potential 
which is needed for the force characteristics. Assuming 
the ideal plastic deformation in the experiments 0, , 
the volume of deformed sample will be constant. All 
considered samples have a same volume, V — 4 nm , 
and the number of electrons N — 236. The linear sizes 
of a sample vary in the range: < L < 13ao for 

slabs, and Lo/10 > a > \f^r s for wires. Here, ag — 
H 2 /me 2 is the Bohr radius. The calculated energy and 
force characteristics corresponding to a slab and wire, 
respectively, are displayed in the left and right parts of 
Figs.l and 2. The wire length, L = Lo + AL, has been 
increased by about seven times. Figure 1(a) shows the 
density of states, D(Ep), of an isolated sample (curve 
1) and the one contacted with reservoirs (curve 2) versus 
the length of wire. For a better demonstration the curves 
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FIG. 1: Characteristics of the elongated sample as functions 
of the length (in units of L/Xfo, where Xfo is the Fermi wave 
length of the reservoir, Xfo = 2n/kFo)- The left side of the 
figure corresponds to a slab while the right one to the wire, 
(a) Density of states at the Fermi level of a isolated sample 
(curves 1) and sample in a contact (curves 2). (b) The work 
function, ionization potential, and contact potential difference 
of the isolated slab and wire. For comparison, the values of 
W (o, •) and IP (□, ■) for rigid (upper points) and self- 
compressed (lower points) spherical cluster of volume 4 nm 3 
[§3 are placed " between " the slabs and wires, (c) The size- 
dependent part of the number of electrons spilling out from 
the sample in a contact with reservoirs, (d) Total kinetic 
energy K of the electrons (left axis) and elongation force F 
(right axis) of the isolated sample. 



2 are slightly shifted up. The peaks of the D(Ep) for a 
wire with square cross-section are more intensive for the 
additionally degenerate subband, k Uy = k nz . 

Figure 1(b) displays the behavior of the electron work 
function and the ionization potential of the isolated sam- 
ples of varying size. The potential difference 8(j> that 
arises at the contact is also indicated. The inequality 
W(L) < Wo is observed to be obeyed over the whole 
range of the considered lengths of the slab and wire. In 
the shortest wire the five subbands appear and four of 
them disappear at elongation. The amplitudes of oscilla- 
tions of the Fermi level are of one or two tenths of an eV. 
This is substantially less than calculated self-consistently 
for the extended thin Al slabs [2^| and wires [23, l3lf . 
As can be seen, there are the ranges of widths, where 



IP < W and IP > W Q . The fact that IP < W 
is rather unexpected. Judging from the empirical fact 
that the work function Wo of the alkali metal is ap- 
proximately equal to one half of the IP of the atom, 
one would expect that the IP of a small solid (inde- 
pendently of the shape of its surface) falls in the range 
W < IP(cluster) < ZP(atom). 

The size correction in W(L) competes with the e 2 /2C 
term in expression 111 Oil for the ionization potential. The 
magnitude of IP(L) depends on the shape of a sample: 
IP(L) > W for L < 43A F0 and IP(L) < W for L > 
43Afo- Note that the positions of the local minima in 
the IP(L) and W(L) curves correspond to the positions 
of the peaks in the density of states of an isolated sample 
(Fig. 1(a)). In the case a = L, we deal with a metal cube 
whose ionization potential must have a value similar to 
that of a sphere. It should be noted that the ionization 
potential depends only on the geometry of the sample 
and is independent of the direction of electron emission 

EH El 

The CPD leads to a noticeable negative shift of the 
potential well depth (in utmost point by about 0.5 eV) 
and this in turn leads (Fig. 1(a)) to a shift in the density 
of states (curve 2) to a zone of the bigger cross section. 

Figure 1(c) shows the size dependence of the number 
of electrons 6N that have left the sample. An increase 
of the cross-section or a decrease in the wire length leads 
to the breaks appearing in the monotonic components of 
S(f>(L) and 8N(L) which correspond to a spherical sam- 
ple. For a of the size close to the atomic diameter the 
SN makes up 20% of the initial number of electrons. The 
account for the contact potential leads to a dependence 
of 5N on the hierarchy of energy levels in an isolated 
sample. This aligns specific breaks in the SN(L) curve 
which the positions of the peaks in the density of states 
D(Ep) of both isolated and contacted samples. The re- 
sults show that some part of the electrons is spilled out 
from the slab (wire) in a contact and the contact can be 
thought of as a sort of the "electron pump" which pulls 
out the electron liquid from or draws it into reservoirs. 
The dipole layers formed in a vicinity of both contacts 
must stimulate additional longitudinal deformation of the 
wire and contribute to a change in its shape and in the 
electron density. The above results allow to calculate size 
dependence of the effective capacitance C c = eSN /5<f> for 
a sample in contact. 

The kinetic energy of electrons and the elongation force 
of the isolated sample are displayed in Fig. 1(d). Their 
counterparts for the sample after junction with reservoirs 
as well as the potential f2 are shown in Fig. 2. Total elec- 
tron kinetic energy K(L) and Ki(L) of metal cube is 
similar to that for semi-infinite metal, (3/5)NEpo- The 
amplitudes of the force oscillations F(L) of the isolated 
slab are several times bigger than those measured for the 
nanowires Q, 0. It is stipulated by a large value of 
the cross-section a x a. Nevertheless, dividing the first 
amplitude of the force by the number of atoms in the 
slab one gets 0.4 nN. This value is smaller than the force 
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FIG. 2: The total kinetic energy Ki, the grand potential Q, 
and the elastic force F\ of elongated slab and wire in a contact. 



acting on a single atom (1.6 nN) in the point contact 
experiments 0, [HI ■ That can be explained by the differ- 
ent dimensionality of a slab and wire. For the reduced 
dimensionality the amplitudes of force oscillations range 
from 1 to 3 nN, which is similar to the measured values. 
The character of the dependence F(L) changes signifi- 
cantly at the wire contact because of depletion of the 
electron subsystem and due to a shift of the spectrum. 
It should be mentioned that in our approach the size de- 
pendencies K (L) and K±(L) are qualitatively different 
from those obtained for an infinitely deep potential well 
00. El The second term appearing in iTH^I only partly 
smoothes out the oscillations in Q(L). This is in contrast 
to the infinitely deep potential well where the oscilla- 
tions in f2(L) disappear. As a result, the amplitudes of 
force F\(L) and their shape change significantly. In our 
more advanced model, the attempt to take into account 
the contact potential difference leads to the appearance 
of spikes in the elastic force versus length which result 
from the fact that for the isolated and contacted wires 
the densities of states are not in line. It is possible that 
in experiment these peaks could not be recognized in the 
background of the thermal fluctuations of the wire-form 
which occur during deformation. With the exception of 
the interval 54 < L/X F0 < 62 (Fig. 2(b)) the calculated 
amplitudes of force agree with the experimental values. 
The difference may be stipulated by not self-consistent 
determination of the CPD for the wires of diameter that 
is comparable with the atomic dimensions and by the 
assumption that L varies continuously. Another reason 
might be that the elastic force is sensitive to the applied 
voltage as suggested in Refs. [33ll34^ . 

The choice of the simplified form of the wire and a 
non self-consistent treatment of the size oscillations in 



the potential U are the disadvantages of the model. The 
free-electron model is not able to describe in detail the 
effect of the atomic rearran gem ents on the magnitude of 
force in the atomic contact [10|. On the other hand the 
obtained qualitative results are quite general and should 
not depend on the symmetry of the problem. It should 
be also noted that for the stabilized jellium model the 
bulk modulus|22ll35j, work function, surface energy, and 
surface stress (^H^ for gold can be described only when 
its valency is assumed to be three. The transition from 
monovalent to trivalent gold can lead to increased, by sev- 
eral times, amplitude of oscillations of the elastic force, 
in analogy with that observed for the trivalent Al [l4| . 
These problems can be connected with an unusual size 
dependence of IP(R) that is measured for spherical gold 
particles: for the very large as well as for the atomic- 
size particles IP > Wo, whereas IP < Wo for parti- 
cles of the intermediate R [37l l3S| |. In our view, this 
could be explained by a reduced valency of some met- 
als which is connected with the decreasing of the size of 
a sample (it reminds the metal-nonmetal transition). A 
similar transition was observed [39] for mercury clusters 
and interpreted in Ref . [4(| . The behavior of the ioniza- 
tion potential measured for aluminium clusters (compare 
Fig. 28(a) in [41|) also points to this effect. Therefore, the 
assumption of monovalcncy of Au in the low-dimensional 
structures is quite reasonable. 

Finally we would like point to the connection between 
our results and the famous phenomenon of the quanti- 
zation of shear modulus for polycrystalline samples that 
was experimentally observed for 57 elements at the room 
temperature |4^. Possibly this effect may be explained 
by a quantization of the force characteristics of the 2D 
electron liquids in a inter-crystallite space of the poly- 
crystal. 

The results of this paper can be summarized as follows. 
The evolution of the size dependence of the ionization 
potential during transformation of the shape of a sam- 
ple from a slab toward a wire is investigated. The size 
dependence of the CPD is calculated in a simple man- 
ner. It is shown that the infinite potential well model of 
a slab |S[ overestimates the value of the CPD by many 
times. For the first time, the effect of the CPD on the 
size oscillations of force was demonstrated. The magni- 
tude of the calculated CPD is significant for the wires of 
subatomic diameters only. An analytical expression for 
the dependence of the potential well depth on the Fermi 
energy for the semi-infinite metallic system is tested for 
low-dimensional structures. [Remarks: Notwithstanding 
that expressions derived in 0, Il6l UtI ITsj successfully 
describe character of the behavior of the electron work 
function (of positive value) for semi-infinite metals and 
ionization potentials for metallic clusters, in our view, a 
usage of the classical electrostatics concept (image force) 
to define the quantum characteristic - work function - is 
not whole appropriate because it cannot unequivocally 
determine the sign of the emitted particle. For example, 
the positron work function is positive for sodium and 
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negative for aluminum in spite of identical action of the supported by the Mianowski Fund (Poland) and partially 
image force on this particle by the Ministry of Education and Science of Ukraine. 
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